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Abstract - In this paper, we explore the use of the Dezert- 
Smarandache Theoiy (DSmT) for seismic and acoustic 
sensor fusion. The seismic/acoustic data is noisy which 
leads to classification errors and conflicts in declarations. 
DSmT affords the redistribution of masses when there is a 
conflict. The goal of this paper is to present an application 
and comparison on DSmT with other classifier methods to 
include the support vector machine(SVM) and Dempster- 
Shafer (DS) methods. The work is based on two key 
references (1) Marco Duarte with the initial SVM 
classifier application of the seismic and acoustic sensor 
data and (2) Arnaud Martin in Vol. 3 with the 
Proportional Conflict Redistribution Rule 5/6 
(PCR5/PCR6) developments. By using the developments 
of Duarte and Martin, we were able to explore the various 
aspects of DSmT in an unattended ground sensor scenario. 
Using the receiver operator cun>e (ROC), we compare the 
methods for individual classification as well as a measure 
of overall classification using the area under the curve 
(AUC). Conclusions of the work show that the DSmT 
results with a maximum forced choice are comparable to 
the SVM. 

Keywords: Information Fusion, DSMT, PCR5, PCR6, 
Area Under the Curve (AUC), SVM, information fusion, 
target tracking, evidential reasoning. 

1 Introduction 

The goal of this paper is to present an application and 
comparison on DSmT with other classifier methods. The 
work is based on two key references (1) Marco Duarte 
with the initial classifier application of the seismic and 
acoustic sensor data [1] and (2) Arnaud Martin in Vol. 3 
for the implementation of the DSmT methods. [2] By 
using the developments of Duarte and Martin, we were 
able to explore the various aspects of DSmT in an 
unattended ground sensor scenario. In the exploration of 
information fusion metrics for classification, there is a 
need to develop metrics of effectiveness that support the 
user’s utility needs [3] and can vary over the sensor types, 
environmental conditions, targets of interest, situational 
context, and users [4]. 

DSmT is an extension to the Dempster-Shafer method of 
evidential reasoning which has been detailed in numerous 
papers and texts: Advances and applications of DSmT for 
information fusion (Collected works), Vols. 1-3 [5], In 2002, 
Dezert [6] introduced the methods for the reasoning and in 
2003, presented the hyper power-set notation for DSmT 
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[7], Recent applications include the DSmT Proportional 
Conflict Redistribution rule 5 (PCR5) applied to target 
tracking [8]. 

The key contributions of DSmT are the redistributions 
of masses such that no refinement of the frame 0 is 
possible unless a series of constraints are known. For 
example, Shafer’s model [9] is the most constrained DSm 
hybrid model in DSmT. Since Shafer’s model, authors 
have continued to refine the method to more precisely 
address the combination of conflicting beliefs [10, 11, 12] 
and generalization of the combination rules [13, 14]. An 
adaptive combination rule [15] and rules for quantitative 
and qualitative combinations [16] have been proposed. 
Recent examples for sensor applications include electronic 
support measures, [17, 18] and physiological monitoring 
sensors [19]. One application of DSmT that has not been 
fully explored is in seismic, magnetic, and acoustic 
classification fusion of moving targets. Kadambe 
conducted an information theory approach [20] and used 
DSmT as integrity constraints [21], but did not take 
advantage of the conflict redistribution. 

Detecting moving vehicles in an urban area [22] is an 
example where DSmT conflicting mass redistribution 
could be helpful [8]. Detecting traffic can be completed by 
fixed ground cameras or on dynamic unattended ground 
vehicles (UGVs). If the sensors are on UGVS, path 
planning is needed to route the UGVs to observe the 
traffic [23, 24] and cooperation among UGVs is 
necessary[25]. The DARPA Grand Challenge featured 
sensors on mobile UGVs observing the environment [26]. 
Mobile sensing can be used to orient [27] or conduct 
simultaneous location and mapping (SLAM) [28]. 

Deployed ground sensors can observe the vehicles; 
however they are subject to noisy measurements from 
environmental conditions. One interesting question is how 
to deploy the fixed sensors that optimize the performance 
of a system. Issues for distributed wireless networks 
(WSNs) include distributed processing, communications, 
and data fusion [29]. In a dynamic scenario, resource 
coordination [30] is needed for both context assessment, 
but also the ability to be aware of impending situational 
threats [31, 32]. For distributed sensing systems, to 
combine sensors, data, and user analysis requires 
pragmatic approaches to metrics [33, 34, 35, 36]. For 
example, Zahedi [37] develops a QOI architecture for 
comparison of centralized versus distributed sensor 
network deployment planning. 

Information fusion has been interested in the problems 
of databases for target trafficability (i.e. terrain 
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information) [38], sensor management [39], and 
processing algorithms [40] from which to assess objects in 
the environment. Various techniques have incorporated 
grouping object movements [41], road information [42, 
43], and updating the object states based on environmental 
constraints [44], Detecting, classifying, identifying and 
tracking objects [45] has been important for a variety of 
sensors, including 2D visual, radar [46], and hyperspectral 
[47] data; however newer methods are of interest for 
ground sensors with ID signals. 

Seismic data provides passive sensing of ground 
vibrations which can be used for motion tracking. Passive 
magnetic sensing can detect hidden objects that might 
indicate intent. Finally, acoustic data can be used for 
signature detection from vehicle engines. [48] The 
DARPA SENSIT program investigated deploying a 
distributed set of wireless sensors along a road to classify 
vehicles as shown in Figure 1. 





Figure 2. Deployed Sensors. [53] 



Track management situational awareness tools receive 
input from sensor feeds (examples include electro-optical, 
radar, electronic support measures (ESMs), and sonar) and 
display this information to a user. User inputs include: 
creation of new objects, such as tracks, contacts and 
targets. Methods to reduce data-to-decisions include: 
fusing multiple tracks into a single track, incorporating 
alerting mechanisms, or visualizing track data common 
operational picture (COP). Sensor and track data can grow 
rapidly as the user desires to keep historical data. 

Our goal is to utilize the DSmT method for the fusion of 
information from seismic and acoustic data in which each 
sensor/classifier is in direct conflict with the other sensor. 
We address (1) intelligent use of the data based on value 
for classification, (2) DSmT sensor data fusion for 
detection, classification, and positional location, and (3) 
metrics to support the sensor and data management as 
supporting a user control. 



Figure 1. SENSIT Data [1] 



2 Location / Detection 



The sensors include acoustic and seismic signals. Given 
the deployed set of sensors, feature vectors were used to 
classify signals based on the data from the seismic and 
acoustic signals. [1] Various approaches include 
combining the data with decision fusion [49], value fusion 
[50], and simultaneous track and identification methods 
[51, 52]. Information theoretical approaches including the 
KL method were applied to the data for sensor 
management [53] as shown in Figure 2. Processing sensor 
data for target classification using acoustic [54, 55] and 
seismic [56] results have been explored in support of 
information theoretical sensor placement [57]. 

Much work has been completed using imaging sensors 
and radar sensors for observing and tracking targets. Video 
sensors are limited in power and subject to day/night 
conditions. Likewise, radar line-of site precludes them 
from observing in the same plane. Together, both imaging 
and radar sensors do not have the advantage of UGSs 
which can power on and off, can work for a long time on 
battery power, and can be deployed to remote areas. 



We desire to track and identify the targets based on the 
sensor reports. In this study, we concentrate on the 
classification of targets which can be used with the 
kinematic/position information for target identification. 

2.1 Sensor Information Management 

The goal is to utilize the UGSs sensors which may be 
acoustic, magnetic, seismic, and PIRoelectric (passive 
infrared for motion detection. With a variety of sensors, 
information fusion can (a) utilize the most appropriate 
sensor at the correct time, (b) combine information from 
both sensors on a single platform, (c) combine results from 
multiple platforms, and (d) cue other sensors in a hand-off 
fashion to effectively monitor the area. Sensor exploitation 
requires an analysis of feature generation, extraction, and 
selection or (construction, transformation, selection, and 
evaluation). To provide track and ID results, we develop 
method or target classification. 

2.2 Sensor Classification 

Sensor exploitation includes detection, recognition, 
classification, identification and characterization of some 
object. Individual classifiers can be deployed at each level 
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to robustly determine the object information. Popular 
methods include voting, neural networks, fuzzy logic, 
neuro-dynamic programming, support vector machines, 
Bayesian and Dempster-Shafer methods. One way to 
ensure the accurate assessment is to look at a combination 
of classifiers. Combination of classifiers [58], could 
include different sensors with classifiers, different 
methods over a single or multiple sensors, and various 
hierarchies of coordinating the classifiers such as Bayes 
nets and Bayesian processing [59], 

Issues in classifier combination methods need to be 
compared as related to decisions, feature sets, and user 
involvement. Selecting the optimal feature set is based on 
the situation and environmental context of which the 
sensors are deployed. An important question for sensor 
and data management is measures of effectiveness. For 
instance, what is the quantification of fusion/decision gain 
using a set of classification methods and placement 
methods? There is a need for a robust combination rule 
that includes the location and detection of the sensors 
subject to the target and environmental constraints. 
Typically, a mobile sensor needs to optimize its route and 
can be subject to interactive effects of pursuers and 
evaders with other targets [60] as well as active jamming 
of the signal [61]. 

Detecting targets from seismic and acoustic data in a 
distributed net centric fashion requires pragmatic 
approaches to sensor and data management. [62] To 
robustly track and ID a target requires both the structured 
data from the kinematic movements as well as the 
unstructured data for the feature analysis. [63] 
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where f] is the interesting and all denominators in the 

equation above are different from zero. If a denominator is 
zero, that fraction is discarded. Additional properties of 
PCR5 can be found in [64], Extension of PCR5 for 
combining qualitative bba’s can be found in [5], Vol. 2 & 
3. All propositions/sets are in a canonical form. A variant 
of PCR5, called PCR6 has been proposed by Martin and 
Osswald in [5], Vol. 2, for combining more than 2 sources. 
PCR6 coincides with PCR5 when one combines two 
sources. The difference between PCR5 and PCR6 lies in 
the way the proportional conflict redistribution is done as 
soon as three or more sources are involved in the fusion. 
For example, let’s consider three sources with bba’s mf.), 
m i(.), and m 3 (.), A D B = 0 for the model of the frame 0, 
and rn , (A ) = 0.6, m 2 (B) = 0.3, and mfB) = 0.1. With PCR5 
the partial conflicting mass mfA) m 2 (B) mfB) = 
(0.6)(0.3)(0.1) = 0.018 is redistributed back to A and B 
only with respect to the following proportions 
respectively: x A PCR5 = 0.01714 and x B PCR5 = 0.00086 
because the proportionalization is [8]: 

Xa PCR5 a- r PCR5 _ mfA)mfB)mfB) 
mi(A ) m 2 (B) mfB) mfA) + m 2 (B) mfE) 



3 DSmT 

Here we use the Proportional Conflict Redistribution rule 
no. 5 (PCR5) and no. 6 (PCR6) and the Dezert- 
Smarandache Probability (DSmP) selections which are 
discussed below. We replace Smets’ rule [10] by the more 
effective PCR5 or eventually the more simple PCR6 and 
replace the pignistic transformation by the more effective 
DSmP transformation to estimate target classification 
probabilities. All details, justifications with examples on 
PCR5 and PCR6 fusion rules and DSmP transformation 
can be found freely from the web in the DSmT compiled 
texts [5], Vols. 2 & 3.. 

3.1 PCR5 and PCR6 fusion rules 

In DSmT (Dezert-Smarandache Theory) framework, the 
PCR5 is used generally to combine the basic belief 
assignment (bba)’s. PCR5 transfers the conflicting mass 
only to the elements involved in the conflict and 
proportionally to their individual masses, so that the 
specificity of the information is entirely preserved in this 
fusion process. Let mf.) and m 2 (.) be two independent 
bba’s, then the PCR5 rule is defined as follows (see [5], 
Vol. 2 for full justification and examples): m PCR5 (0) = 0 
and VW £ 2 s \ {0}, where 0 is the null set and 2 s is the 
power set: 



x A PCR5 _ x r PCR5 _ 0.018 

thatls 0.6 (0.3)(0.1) 0.6 + 0.03 ~ °- 02857 

thus x A PCR5 = 0.60 (0.02857) = 0.01714 
x B PCR5 = 0.03 (0.02857) = 0.00086 



With the PCR6 fusion rule, the partial conflicting mass 
mfA) m 2 (B) mfB) = (0.6)(0.3)(0.1) = 0.018 is 
redistributed back to A and B only with respect to the 
following proportions respectively: x A PCR6 = 0.0108 and 
x b PCR6 = 0.0072 because the PCR6 proportionalization is 
done as follows: 



x A 



*B:2 



*B:3 



m,(A) m 2 (B ) mfB) 



»U(A ) m 2 (B ) mfB) 
m\(A ) + m 2 (B) + mfB) 



that is 



x A 



0.6 



XB;2 



Xb;3 



0.3 



0.1 



0.018 

0.6 + 0.3 + 0.1 



0.018 



thus 



x A =(0.6) (0.018) = 0.0108 
x B / CR6 = (0.3) (0.018) = 0.0054 
x b ” 3 pcr6 = (0.1) (0.018) = 0.0018 



81 



and therefore with PCR6, one gets finally the following 
redistributions to A and B: 

jc a pcr6 = (0.6) (0.018) = 0.0108 



PCR6 

Xg 



~ ^B,2 



PCR6 



+ 



x B / CR6 = 0.0054 + 0.0018 = 0.0072 



derived by choosing e as a very small positive number, say 
e = 1/1000 for example in order to be as close as we want 
to the highest value of the PIC. The smaller e, the 
better/bigger PIC value one gets. When e = 1 and when the 
masses of all elements Z having C(Z) = 1 are zero, 
DSmP e =i(.) = BetP(.). 



From the implementation point of view, PCR6 is simpler 
to implement than PCR5. For convenience, Matlab codes 
of PCR5 and PCR6 fusion rules can be found in [5]. 

3.2 The DSmP Transformation 

The DSmP probabilistic transformation is an alternative to 
the classical pignistic transformation which allows us to 
increase the probabilistic information content (PIC), i.e. to 
minimize the Shannon entropy, of the approximated 
subjective probability measure drawn from any bba. 
Justification and comparisons of DSmP(.) with respect to 
BetP(.) and to other transformations can be found in 
details in [65, 5 Vol. 3, Chap. 3]. 

BetP: The pignistic transformation probability, denoted 
BetP, offers a compromise between maximum of 
credibility Bel and maximum of plausibility PI for 
decision support. The BetP transformation is defined by 
BetP(0) = 0 and VXe G 0 \ {0}by 

BetPI. X) = V MXgPmp) W 

Y£G* V ' 

where G 0 corresponds to the hyper -power set including all 
the integrity constraints of the model (if any ). G® = 2® if 
one adopts Shafer’s model for 0 and G® = D® 
(Dedekind’s lattice) if one adopts the free DSm model for 
0 [5], Cjy[(F) denotes the DSm cardinal of the set Y, 
which is the number of parts of Y in the Venn diagram of 
the model M of the frame 0 under consideration [5, Book 
1, Chap. 7], The BetP reduces to the Transferable Belief 
Model (TBM) when G® reduces to classical power set 2® 
when one adopts Shafer’s model. 

DSmP transformation is defined by DSmPe(0) = 0 and 
VX e G® \ {0} by: 



DSmP t (X) 



Y m(Z)+e-C(X nY) 

ZCXnY 

C(Z)=1 

V ^ m(Y) 

yTZ* Y rn(Z) + e-C(Y) 

ZCY 
C(Z )= 1 



(3) 



where C(X D Y) and C(y) denote the cardinals of the sets 
If! 7 and Y respectively; e > 0 is a small number which 
allows to reach a highest PIC value of the approximation 
of m{.) into a subjective probability measure, and Z is the 
new evidence. Usually e = 0, but in some particular 
degenerate cases, when the DSmP e =o(.) values cannot be 
derived, the DSmP E >o values can however always be 



4 Seismic/ Acoustic Example 

We use the SEN SIT data which was described previously. 
To perform the data analysis we use data mining [66] 
techniques such as a support vector machine (SVM) [67, 
68] to process the unstructured data. Through analysis, we 
can determine the optimum use of the data given 
environmental conditions (i.e. obscurations) and sensor’s 
capabilities to detect a moving target. 

Figure 3 shows the methodology of comparison. A key 
comparison is made between combining all the acoustic 
and seismic data together for testing and training via the 
SVM versus using the outputs from the acoustic and 
seismic data separately from which conflicts in 
classification are detected and sent to DS and DSmT 
processing. 




For the data, we have 50 seismic and 50 acoustic 
measurements at each time step for 10,000 time steps. We 
trained the SVM over the first 500 measurements and ran 
the study over a different 500 time steps. With limited 
explanation of the all the data, we note that contextual 
factors could vary the measurement quality over the 
different time intervals we chose to test and evaluate, but 
we are unable to verify differences over the intervals. The 
combined decision m(») is from two sources with three 
target choices m x ( A), m x ( B), m x (C), where X is acoustic, 
seismic, or combined decisions. We implement as a 
probability vector (e.g. [/«a(1), «a(2), «a(3)], for acoustic 
with target labels of 1, 2, 3). We then use the DSmP with a 
cardinality of 3. As per (l)-(3), the combined results use 
[m A (l), m A (2), m a (3) m s ( 1), m s ( 2), m s ( 3)], where source A 
is acoustic and source S is for seismic and the outputs are 
Xi PCR6 , X2 PCR6 , and x 3 PCR6 . It is noted that we have not used 
the raw measurements, but the aggregate decision trained 
from the SVM over the 50 measurements. 

Using the results of the analysis, we compared the 
decisions versus the truth data over all decisions for 100 
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ROC curve 



trails and plotted a receiver operator curve (ROC). The 
ROC plots the probability of detection ( P D ) or true 
positive rate (TPR) versus probability of false alarm (P F a) 
or False Positive Rate (FPR). To further gather a 
meaningful measure we plot the area under the curve 
(AUC) to compare the different cases. 

4.1 Data Processing 

We compare two cases of (1) processing the data 
separately and (2) jointly processing the acoustic and 
seismic results Figure 4 shows the case of the acoustic 
results where the AUC for each target is calculated as an 
overall performance criterion. 



ROC 




Figure 4. Acoustic Results. 



Figure 5 plots the seismic results. Note that for the data 
set, the seismic results have a lower probability of false 
alarms for target 3 and target 2; however, target 2 exhibits 
more confusion. 



ROC curve - Seismic 




Figure 5. Seismic Results. 



Next we explore the case of the joint seismic and acoustic 
data management and utilize SVM for classification, 
shown in Figure 6. Note the false alarm reduction (FPR) 
which is desired by users. 




Figure 6. Combined Results 



In general, the joint analysis supports better decision 
making with improved confidence (TPR or P D ) for a 
constant false alarm rate, target location accuracy from 
joint spatial measurements, and timeliness (i.e. target ID 
decision made with two modalities required fewer 
measurements than that of a single modality). 

4.2 Application of DS 

In this section, we show the results of the application of 
DS methods. Given a training and prediction results in a 
combined probability, we have for target 1, target2, and 
target3 a vector of m = [ni\ m 2 m 2 \. Based on the 
prediction results from the SVM, there are many conflicts 
of the sensor decision based on the maximum probability. 
When a conflict occurs, it would be better suited to 
acknowledge the conflict and then redistribute the 
probabilities based on a set notation. In this case, the focal 
elements are 0 = [fl, ..., (9 7 ] = [‘l’, ‘1(12’, ‘2’, ‘2fj3\ ‘3’, 

‘ lfl3 ’, ‘1(12(13’]. Using the analysis by Martin, we 
conduct an analysis over the set criterion. Figure 7 shows 
that a significant reduction false alarms (at a low FPR the 
TPR remains high). 



ROC curve - DS 




Figure 7. DS allowing for set declarations. 



The shift of the curve to the left (versus Figure 6 of the 
SVM) demonstrates that DS is better at low _P F A . Flowever, 
if we compare the AUC for all results, then the SVM 
performs better. Thus, there is a trade off when for using 
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ROC curM3 - DSMT1 



DS for increasing the P» at low P F A versus the SVM for 
overall classification analysis. 

To explore a comparison of approaches, we utilized the 
bba and forced the evidential reasoned to choose a single 
target. From this analysis, the AUC improves in 
comparison to the SVM approaches which are a forced 
choice analysis. Figure 8 plots the DS (for one target 
designation). 



ROC curve - DS1 




Figure 8. DS Single Target Detection. 



4.3 Application of DSmT 

DSmT, as described above, improves on the methods of 
DS conflict redistribution. In this example, the bba were 
computed from the SVM and calculated with (l)-(3). 
DSmT classification was comparable to the complete 
SVM fusion analysis for the AUC, with improvements at 
the low _P F a- 

Figure 9 presents the DSmT results for set declaration 
and Figure 10 shows the case of a forced target choice 
from the DSmT. 



ROC curve - DSMT 




Figure 9. DSmT allowing for set declarations. 



From Figure 9, we can see that the set-based approach 
improves the detection for low false alarm rates. At high 
false alarm rates or FPR (0.3 — 0.9), the detection 
probability for DSmT is lower than the SVM results. 
Using the maximum of the target bba, shown in Figure 10, 
provides an analysis threshold that renders the DSmT 
comparable to a SVM. 




Figure 10. DSmT Single Target Detection. 



In the table below, we look at the entire analysis using the 
area under the ROC (AUC) as a key metric in the analysis. 
There are cases in which the maximum AUC and 
minimum AUC are improved but the overall analysis 
(Total AUC) varies. We see from the comparison three 
results : (1) that the DSmT improves upon the DS 
methods, (2) information fusion improves the results, and 
(3) the DSmT without redistribution of the conflict for 
each measurement can obtain results similar to the SVM. 



Table 1: AUC Comparisons of SVM, DS, and DSmT 



Method 


Min AUC 


Max AUC 


Total AUC 


A-SVM 


0.786 


0.821 


2.401 


S-SVM 


0.696 


0.844 


2.335 


C-SVM 


0.791 


0.851 


2.472 


DS 


0.671 


0.742 


2.141 


DS1 


0.738 


0.833 


2.371 


DSMT 


0.728 


0.751 


2.224 


DSMT1 


0.760 


0.855 


2.440 



A — Acoustic, S-Seismic, C-Combination 



Table 1 has three motivations. In the first look, for a 
single target detection; the SVM alone (run over all the 
data) does perform slightly better than the DSmT 
information fusion case. Second, from the results we 
advocate the use of DSmT over DS. Third, we see that 
while DSmT is more complex, it offers further capability 
to look at each measurement and redistribute the conflict 
versus using all of the data in a force SVM case. 

If we compare Figure 6 (SVM combined) and Figure 10 
(DSmT max choice), we see that all curves in the DSmT 
have better results at a low P F A rate over all targets; 
whereas the SVM has lower performance for target 1 and 
target 2 and a better classification of target 3. The 
aggregate AUC measure provides a global comparison; 
which advocates the SVM; however, DSmT is more 
conservative through the proportional conflict 
redistribution which improves the SVM analysis over the 
low P F A rate results. The P F a is important for decision 
support as operators are interested in the P D at low P FA . 
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5 Conclusions 

We have explored DS and DSmT methods for seismic and 
acoustic information fusion. The goal of the paper was to 
extend the existing techniques presented by Martin and 
Durate for further demonstration of how to deal with 
conflicting information. From the initial results, the use of 
DSmT can be tailored to the seismic and acoustic sensors 
for on-line high conflict mass redistribution of decision 
outputs as the sensors measure different target 
phenomenologies and might not always be globally 
optimized as with the SVM. We utilized a Bayesian basic 
belief assignment (bba) with only singleton as focal 
elements which from the P vectors of the target 
probabilities. Future work will use non-Bayesian 
approaches to get the bba’s, further assessment of the 
conflict redistribution, set-based decision logic, and 
application to evidential-based simultaneous tracking and 
identification [69, 70, 71] Additional applications include 
group tracking [72, 73, 74, 75] and video tracking [76]. 

Information theoretic measures [77] and tracking 
analysis [78] can support the sensor and data management 
as well as determine the Quality of Information and 
Quality of Service needs. Use of the Area Under the Curve 
(AUC) provides decision support for situational awareness 
for command and control from which we can extend to 
higher dimensions [79], Various other sources of soft data 
(human reports) can be combined with the hard (physics- 
based sensing) [80] to update the sensor management, 
situation awareness [81], and reporting of the situation 
based on the context and the needs of users such as 
measures of effectiveness for mission support. 
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